
Figure 2: The distribution of pulls SA/afor a measurement 
of A using the method of moments on samples 0/I6OO events. 
The curve shows the result of a fit to a gaussian distribution, 
with (6A) = -0.005 ± 0.026 and a = 1.012 ± 0.020. 



6 A Monte Carlo evaluation of the method 

With a simple Monte Carlo program, the method of moments was compared to the unbinned 
likelihood method. The probability of the likelihood method was taken to be proportional to 
N + (C cos 9) as given by eq. (2). The distribution as generated in the Monte Carlo could be 
different from this, to simulate an angular dependence in the efficiency or the purity. For different 
numbers of Monte Carlo events one can determine the measured value of A for both methods 
and the estimated standard deviations. Such a Monte Carlo measurement can be repeated many 
times. This makes it possible to determine the real standard deviations as well. 

An internal check of both methods and of the program itself is to look at the distribution of 
the 'pulls': the difference between the input asymmetry and the measured values, divided by 
the estimated standard deviation. As it should be, the distribution of the pulls can be fitted very 
precisely by a gaussian centred at zero and with a width equal to 1 to an accuracy of about 1% 
(see e.g. figure 2). 

It turns out that the method of moments is equivalent to the unbinned likelihood method: on 
average, the standard deviations are equal. Both methods are very robust against uncertainties 
in the detection efficiency. The difference between the two methods is less than a few percent 
of the error estimate for measurements using numbers of events ranging from 50 to 256000. 

7 Conclusion 

The method of moments developed here is equivalent to the unbinned likelihood method, but 
simpler to use and it gives useful analytical formulas for the expected statistical errors in certain 
simple but quite realistic cases. This is as good as one can get, since it is well known that the 
unbinned likelihood method is, nearly by definition, equivalent to the best possible method [1]. 
Both methods can be refined to include the effects of the efficiency and the background. 

References 

[1] One of the many good textbooks on statistical methods is: W.T Eadie, D. Dryard, EE. James, 
M. Roos, B. Sadoulet, ' Statistical Methods in Experimental Physics', North-Holland Pub- 
lishing Company, 1971 
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e min (degrees) 9 min (degrees) 

(a) fixed number of events inside the acceptance (b) fixed number of Z° decays 

Figure 1: The dependence of the standard deviation a on the polar acceptance. Figure (a) 
visualises the coefficient ofl/N in equation ( 1 7). This gives the expected statistical error for a 
given number of detected events inside the angular acceptance. Figure (b) shows the angular 
dependence of a 2 from equation (18), normalised to the case of full acceptance (9 m m = 0). 
This second plot shows the relative importance of different angular regions of a detector. For 
detectors at LEP, the borders between the 'barrel' and 'endcap' regions are near to 40° and 
140°. 



/ -> (1 + D) = f. We can call D the 'dilution', while R = 1/(1 + D) is the 'purity': the 
fraction of the original signal in the total sample. After a renormalization of the distribution one 
can see that the detected asymmetry decreases asA^A/(l + D) = RA. Let the measured 
asymmetry before correcting for the background be given by A ± s/ \/N. The accuracy of the 
corrected result will now depend on the number of events, the purity R and on its uncertainty 
SR: 

A = RAq CT 2 = _f_ + {AoSR) 2 (19) 

This shows that it is more important to have a high and well known purity R 1 than to have 
a high detection efficiency (large N). 

Of course it is also possible that the background is not not proportional to /(cos 9), or is not 
even a symmetrical function of cos 9. In this case the symmetry properties of / (and possibly 
also g) can be modified. Such behaviour of the background can be taken into account by using 
the more complete relation (13). However, this is only possible after a precise determination of 
the background distribution. 
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3 The accuracy and the angular acceptance at the Z° 

Equation (15) provides a simple formula for the expected accuracy of a measurement of A FB 
at the Z° with a reduced angular acceptance, given by | cos(0) | < cos # m in • Since we know that 
A < 0.1, we can now also drop the 0(A 2 ) term in the error estimate. The function h = g/h for 
the distribution (2) is equal to: 

h(cose) = 8 -. ; Z ° S9 (16) 
o (1 + cos^ 9) 

After normalising distribution (2) inside the acceptance, the standard deviation of A is found to 
be: 

a (^\2 ^ 1 _ _3_ COS 6> m in (3 + COS 2 6> m in) 1 

(h 2 ) ■ N tot 64 cos 6> min - arctan cos 6> m in N tot 
Here N tot stands for the number of events inside the angular acceptance. Equation (17) is 
visualised in figure 1(a) and can be used to find the expected statistical error for a given angular 
acceptance and a given number of events inside the acceptance. 

To understand the importance of a large polar acceptance when measuring Afb, one should 
also express the accuracy in terms of a the number of Z° decays N%. In this case distribution (2) 
should be normalised for the full acceptance: 

v(A) 2 ot — • — -— (18) 

cos # m i n — arctan cos # m i n I\z 

This dependence is visualised in figure 1(b), from which it is easy to see that for a measurement 
of Apb the angular acceptance of the endcaps (20° < 9 < 40° and 140° < 9 < 160°) is 
approximately as important as that of the barrel region (40° < 9 < 150°). In other words: 
adding the information of the endcap-detectors is approximately equivalent to doubling the 
statistics. 



4 The experimental detection efficiency 

The symmetry properties of / and g can in principle be altered by asymmetries in the detection 
efficiency. Corrections due to such asymmetries can be assumed to be small, because of 
cancellations in the ratio of particles to anti-particles. The derivation of the moments started 
by only using ratios in which the detection efficiency cancels in every angular bin separately. 
As long as the detection efficiency does not depend on the 'charge' C, the answer will still be 
correct, although it could be that the error estimate is slightly biased by the angular dependence 
of the weights. If the angular dependence of the efficiency E(cos9) can be determined, this 
information can be used to make the corrections / — >■ Ef and g — >■ Eg. Then eq. (15) is again 
valid. If the detection efficiency is very asymmetric in cos 9, it could be necessary to use the 
full equation (13), including the correction oc A 2 . 

5 An experimental background 

When the signal is not pure, the background can often be assumed to be symmetrical in cos 9. 
Often such backgrounds give an extra contribution to the symmetrical part of distribution (3): 
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where the weights w e depend on 8A e , the standard deviation of A e : 



(xa ? = (&e)2 = 4iVg+iVe ~ = l^i = 1 - A " h2 (9) 

{ B) h(9)~ 2 N$h? N e h 2 N e h 2 } 

To simplify the denominators of the weights in expression (8) we substitute the measurement 
Ae by the constant A. Then definitions (6) can be used to rewrite sums over 9 in terms of sums 
over events i. Subsequently, definition (5) shows how to express the resulting sums as averages. 

Neaph C t h^) M / Ch \ 

^ WeAe = ? T^W = 1 - A 2 h(9i) 2 = Ntot \rW/ 00) 

_ N e h 2 h(9 t ) 2 I h 2 \ 

^ We ~ ^l-A 2 h 2 ~^ t l-A 2 h(9 t ) 2 ~ tot \l-A 2 h 2 / { } 

In the weighted mean the factors N tot cancel to give: 

A = (i^h 2 ) ' (r=W) ±N ^ l2 (i-A 2 h 2 ) (12) 

For sufficiently small \A\ we can expand the right-hand side of eq. (12) in orders of A 2 , 
dropping all terms 0(A 4 ): 

a = ^ + A Ah*)(Chyvh)(h*) ± N;1 , 2 {{h2) + A2( , 4)) - 1/2 + 0(A4) (13) 

This equation can be solved by substituting the first order solution A = (Ch)/(h 2 ) on the 
right-hand side. 

It can be shown that the second term on the right-hand side of equation (13) vanishes. This 
follows from identities for the distribution (3) that can be derived after invoking the symmetry 
properties (4) of / and g: For positive integers k and n, we have: 

J 1 f k (cos9)g 2n+1 (cos9)dcos9 = (Ch 2k+1 ) = A(h 2k+2 ) (14) 

Substituting the result for A in the error estimate, equation (13) is now simplified to: 

A = ^±N;: t /2 ((h 2 ) + (Ch) 2 (h')/(h 2 ) 2 )- 1/2 + 0(A 4 ) (15) 

This formula already includes effects of a reduced acceptance, or even the effects of a detection 
efficiency that is non-uniform in cos 9, but symmetric in the charge C. 

It seems as if nothing is left of the binning we started with. However, the above derivation 
did rely on the existence of sufficient statistics in each bin of cos 9 for efficiency effects to cancel 
in the measured ratio a e = (Ng~ — N e ~)/(N e + + Ng~). It is thus advisable to reject values of 9 
with a very low efficiency. 
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weights, and express the result in terms of moments. To do this we will switch from taking the 
average over bins to taking the average over events. 

The mathematical derivation in this section will actually be done for a more general angular 
distribution. Distribution (2) it is a function of the combination C cos 9. The most general 
distribution with this property can be described by (C = ±1): 

N c (cos 9) oc f (cos 9) + AC g (cos 9) (3) 

where the functions / and g have to satisfy the following symmetry properties: 

f(-x) = f(x) g(-x) = -g(x) (4) 

such that g(C cos 9) = Cg (cos 9) and f(C cos 9) = /(cos 9). The asymmetry parameter A is 



x. 



normalised by requiring J_ t \g(x)\dx = J_ t f(x)d. 

The moments of the distribution of a given observable are defined as the averages of powers 
of this observable. In the case at hand an observable can be any function of the charge C and 
the angle 9, taking on the value Oi = 0(0;, C'i) for each event i. Let Oe be the average of 
inside an angular bin centred at 9. Now the average (0) can be defined by a sum over events, 
but also by a sum over the averages Oe : 

(O) = E°(^/^tot = £O e -iViV tot (5) 

all i e 

The subscript 9 denotes that the sums are over events i in the bin centred at the angle 9. Nq is 
the number of events in the bin at 9, while N tot is the total number of events. The n-th algebraic 
moment of is defined to be (O n ), while the corresponding central moment is ([0 — (0)] n ) . 

To derive the optimal moments to use, we will start by examining the data in different bins 
of 9. In every bin of 9, the following quantity can be measured: 

= N+ - Nf = Zi at e Q 

In such a ratio the detection efficiency cancels for each bin separately. 

Equation (3) relates a measurement a e in a bin centred at 9 directly to the asymmetry 
parameter A: 

a e = Ah(9) , with h(9) = g(9)lf(9) (7) 

To derive an expression for A in terms of moments we will determine the weighted average 
of the measurements of A in all the bins. When the theoretical distribution is used to obtain 
this value, the result is an equation defining A implicitly. By making some simple and plausible 
assumptions, one can subsequently express the weighted average and its statistical accuracy in 
terms of moments and get rid of the binning. 

Let A e = a e /h(9) denote the measurement of A in a given bin. The weighted mean A of 
the values A e and its error are defined by: 

A = %^±(E^r 1/2 w e = (6Ae)- 2 (8) 
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1 Introduction 



Parity violation occurs both in the production and decay of the Z° boson at an e + e~ collider. 
This results in an asymmetry in the angular distribution of the produced fermions: the forward- 
backward asymmetry, Afb- Equivalently, one can measure an asymmetry between particles 
and anti-particles: the charge asymmetry, A c . 

The conventional way to define this asymmetry is by counting particles and anti-particles in 
the forward and backward hemispheres as defined with respect to the direction of the incoming 
electron: 

/v + — /v + /v + — n~ 

a forward backward forward forward \ / 1 \ 

Afb ~ — — t^f — - Ac (1) 

forward 1 backward forward 1 forward 



The equality of the two asymmetries can be understood as a consequence of CP invariance. 

Although the definition of these asymmetries is clear, a direct application of this definition 
is not the best way to measure them. The reason is that the data in different angular regions 
should be weighted differently, according to the changing sensitivity. To see this, we have to 
look at the explicit angular distribution. 

A 'charge', C = ±1, can be assigned to the fermion and the anti-fermion of the pair created 
in the Z° decay. The angular distribution can be given as follows: 

N c cx-( l + cos 2 9) + AC cos 9 (2) 

8 

where the polar angle 9 is defined with respect to the incoming electron. This distribution 
depends only on the combination C cos 9. It is easy to check that the asymmetry parameter A is 
equal to the definitions given above for A FB and A c . 

From this distribution it can be seen that the sensitivity to A increases with increasing values 
of | cos 9\. The standard way to optimise measurements is to use an unbinned log(likelihood) 
fit [1]. This is a method that will give the best possible accuracy. In this article we will derive 
a simpler method that is based on the moments of a special variable. This method of moments 
is equivalent to a likelihood fit, but it gives explicit analytical formulas for both the asymmetry 
and its statistical error. The dependence on the angular acceptance of the experiment can thus 
be given analytically. 

We will first derive the method of moments for a more general case. Subsequently, the 
interesting formulas for our explicit problem of A FB = A c will be given. By construction, the 
result of the method will not change due to the experimental detection efficiency, even if the 
efficiency varies (slowly) as a function of the angle 9. It is also possible to treat the effects of a 
background. Finally a Monte Carlo program was used to compare this method to an unbinned 
likelihood fit. 



2 The method of moments 

Expressing the asymmetry and its statistical accuracy directly in terms of moments has the 
advantage of avoiding the requirement to bin the data. However, to derive the best suited 
moment we will start from the fact that the charge asymmetry can be measured in any bin of 9 
separately. Subsequently, we will combine the measurement of many bins in 9 with the proper 
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Abstract 

Parity violation at LEP or SLC can be measured through the charge asymmetry. An 
optimal method of moments is developed here to measure this asymmetry, as well as similar 
asymmetries. This method is equivalent to the likelihood fit. It is simpler in use, as it gives 
analytical formulas for both the asymmetry and its statistical error. These formulas give the 
dependence of the accuracy on the experimenatal angular acceptance explicitly. 



